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ABSTRACT 

Based on the two-parameter family nature of X-ray clusters of galaxies obtained in a separate 
paper, we discuss the formation history of clusters and cosmological parameters of the universe. 
Utilizing the spherical collapse model of cluster formation, and assuming that the cluster X-ray 
core radius is proportional to the virial radius at the time of the cluster collapse, the observed 
relations among the density, radius, and temperature of clusters imply that cluster formation 
occurs in a wide range of redshift. The observed relations favor the low-density universe. 
Moreover, we find that the model of n ~ — 1 is preferable. 

Subject headings: cosmology: theory — clusters: galaxies: general — X-rays: galaxies 

1. Introduction 

Galaxy clusters are the largest virialized objects in the universe and provide useful cosmological probes, 
since several properties of clusters are strongly dependent on the cosmological parameters. For example, 
the statistics of X-ray clusters can serve as an excellent probe of cosmology. The abundance of clusters 
and its redshift evolution can be used to determine the cosmological density parameter, £Iq, and the rms 
amplitude of density fluctuations on the fiducial scale 8/i _1 Mpc, as (e.g. White, Efstathiou, & Frenk 1993 
; Eke, Cole, & Frenk 1996 ; Viana & Liddle 1996 ; Bahcall, Fan, & Cen 1997 ; Fan, Bahcall, & Cen 1997). 
Moreover, temperature and luminosity function of X-ray clusters is also used for a cosmological probe. 
Taking account of the difference between cluster formation redshift and observed redshift, Kitayama & 
Suto (1996) computed a temperature and luminosity function semi-analytically; comparing the predicted 
temperature (T) and luminosity (Lx) function with the observed ones, they conclude that Slo ~ 0.2 — 0.5 
and h ~ 0.7. However, one tenacious problem in such investigations is the discrepancy of Lx — T relation 
between observations and simple theoretical prediction. This relation should also be an important probe of 
the formation history and cosmology. 

In a separate paper (Fujita & Takahara 1999; hereafter Paper I), we have shown that the clusters 
of galaxies populate a planar distribution in the global parameter space (log p$, log R, logT), where po is 
the central gas density, and R is the core radius of clusters of galaxies. This 'fundamental plane' can be 
described as 

X = p A7 R°- 65 T- om = constant . (1) 

We thus find that clusters of galaxies form a two-parameter family. The minor and major axes of the 
distribution are respectively given by 

Y = ^-39^0.46^.80 ; (2 ) 



-2 - 



Z = p o-™ R -o.ei T -o.o39 (3) 

The scatters of observational data in the directions of Y and Z are AlogY" = 0.2 and AlogZ = 0.5, 
respectively (Paper I). The major axis of this 'fundamental band' (Z) is nearly parallel to the logi? — log/? 
plane, and the minor axis (Y) describes the Lx — T relation. 

In this Letter, we discuss cosmological implications of the relations we found in Paper I, paying a 
particular attention to the two-parameter family nature of X-ray clusters and to the difference between 
cluster formation redshift and observed redshift. In §||, we use spherical collapse model to predict the 
formation history of clusters of galaxies, and in §|J we predict the observable distribution of X-ray clusters. 



2. Formation History of Clusters of Galaxies 

In order to explain the observed relations between the density, radius, and temperature of clusters of 
galaxies, we predict them for a flat and an open universe theoretically with the spherical collapse model 
(Tomita 1969; Gunn & Gott 1972). For simplicity, we do not treat vacuum dominated model in this paper. 
For a given initial density contrast, the spherical collapse model predicts the time at which a uniform 
spherical overdense region, which contains mass of M v ; r , gravitationally collapses. Thus, if we specify 
cosmological parameters, we can obtain the collapse or formation redshifts of clusters. Moreover, the model 
predicts the average density of the collapsed region p v j r . 

In Paper I, we showed that the observed fundamental band are described by the two independent 
parameters M v j r and p v i r . In particular, the variation of p v ir is basically identified with the scatter 
of Z. Since the spherical collapse model treats /9 v ir and M v ; r as two independent variables, it can be 
directly compared with the observed fundamental band, as long as we assume that the core radius and 
the mass of core region are respectively a fixed fraction of the virial radius and that of the virial mass 
at the collapse redshift, as is adopted in this paper. Although the model may be too simple to discuss 
cosmological parameters quantitatively, it can plainly distinguish the relations between the density, radius, 
and temperature in a low-density universe from those in a flat universe, as shown below. 

For the spherical model, the virial density of a cluster is A c times the critical density of a universe at 
the redshift when the cluster collapsed {z co \\). It is given by 

A I \ A J-?/ \ 2 a ^OPcrit,o(l + ^coll) 3 ,.s 

Si(^G0llJ 

where is the cosmological density parameter, and E(z) 2 — r2 (1 + z) 3 /Sl(z). The index 

refers to the values at z = 0. Note that the redshift-dependent Hubble constant can be written as 
H(z) = 100hE(z) km s _1 Mpc -1 . We fix h at 0.5. In practice, we use the fitting formula of Bryan & 
Norman (1998) for the virial density: 

A c = 18tt 2 + 60x - 32x 2 , (5) 

where x = £l(z co \\) — 1. 

It is convenient to relate the collapse time in the spherical model with the density contrast calculated 
by the linear theory. We define the critical density contrast 6 C that is the value, extrapolated to the present 
time (t — t ) using linear theory, of the overdensity which collapses at t = t co ii in the exact spherical model. 
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It is given by 



*c(*cou) = ^D(t ) 



( tg 
Vcoll 



2/3' 



(O < 1) 



3(12tt) 2 / 3 / t 



20 



\icoii 



2/3 



(O = 1) 



(6) 



(7) 



(Lacey & Cole 1993), where D(t) is the linear growth factor given by equation (A13) of Lacey & Cole 
(1993) and i n = 7rff " 1 fi (l - fi ) _3/2 - 

For a power-law initial fluctuation spectrum P oc k n , the rms amplitude of the linear mass fluctuations 
in a sphere containing an average mass M at a given time is 5 oc M~( n+3 ^ 6 . Thus, the virial mass of 
clusters which collapse at i co ii is related to that at to as 

, -I -6/(n+3) 



M vir (i co ii) = M vir (t ) 



Sc(tcou) 



Sc(t ) 



(8) 



Here, M v - 1T (to) is regarded as a variable because actual amplitude of initial fluctuations has a distribution. 
We relate t = t co u to the collapse or formation redshift z co u, which depends on cosmological parameters. 
Thus, M v i r is a function of z co n as well as M v i r (to). This means that for a given mass scale M v ; r , the 
amplitude takes a range of value, and thus spheres containing a mass of M v ; r collapse at a range of redshift. 
In the following, the slope of the spectrum is fixed at n — — 1, unless otherwise mentioned. It is typical in 
the scenario of standard cold dark matter for a cluster mass range. 



The virial radius and temperature of a cluster are then calculated by 



/3M V 



1/3 



T ■ 

1 vi r 



fima GM V 



(9) 



(10) 



3fc B r vir 

where /x(= 0.6) is the mean molecular weight, mjj is the hydrogen mass, fce is the Boltzmann constant, and 
G is the gravitational constant. 



3. Results and Discussion 

Since equations (|4|) and (^) show that p vlY and M v i r are the functions of z co \\ for a given M V i r (to), 
the virial radius r v ; r and temperature T v ; r are also the functions of z co \\ for a given M v { T (to) (equations 
H and [fTof ). Thus, by eliminating z co \\, the relations among them can be obtained. Since observational 
values reflect mainly the structures of core region while the theory predicts average values within the 
virialized region, we must specify the relation between the observed values (po , R, T) and theoretically 
predicted values (p v irj -Rvir, T v - a ). Since we assume that mass distribution of clusters is similar, r v - lr oc R and 
T v i r = T, emphasizing that r V j r is the virial radius when the cluster collapsed (see Salvador-Sole, Solanes, & 
Manrique 1998). In this case, the typical gas density of clusters has the relation po oc /p v ; r , where / is the 
baryon fraction of the cluster. Since it is difficult to predict / theoretically, we assume / oc M°; r 4 , which is 
consistent with the observations and corresponds to X ~ constant (Paper I). For definiteness, we choose 
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Figurc 1 shows the predicted relations between fp v iv, ?~vir> and T v ; r for f^o = 0.2 and 1. Since we are 
interested only in the slope and extent of the relations, we do not specify M v ; r (to) exactly. Moreover, since 
M V i T (to) has a distribution, we calculate for M v i r (i ) = 10 16 Mq and 5 x 10 14 Mg. The lines in Figure 1 
correspond to the major axis of the fundamental band or Z because they are the one-parameter family of 
p v i r . The width of the distribution of M v i r (to) represents the width of the band or Y. The observational 
data of clusters are expected to lie along these lines according to their formation period. However, when 
Oo = 1, most of the observed clusters collapsed at z ~ because clusters continue growing even at z = 
(Peebles 1980). Thus, the cluster data are expected to be distributed along the part of the lines close to 
the point of z co \\ — (segment ab). In fact, Monte Carlo simulation done by Lacey & Cole (1993) shows 
that if f2o = 1, most of present clusters (M v ; r ~ 10 15 Mq) should have formed in the range of z co a < 0.5 
(parallelogram abed). When Q — 0.2, the growing rate of clusters decreases and cluster formation gradually 
ceases at z <^ 1/fio — 1 (Peebles 1980). Thus, cluster data are expected to be distributed between the point 
of z co n = and z co \\ — l/^o — 1 an d should have a two-dimensional distribution (parallelogram ABCD). 

The observational data are also plotted in Figure 1. The data are the same as those used in Paper I. 
For defmiteness, we choose r vir = 8R and /p v ir = 0.06po- The figure shows that the model of O = 0.2 is 
generally consistent with the observations. The slopes of lines both in the model of f^o = 0.2 and flo = 1 
seem to be consistent with the data, although the model of ilo = 0.2 is preferable. However, the model of 
flo = 1 is in conflict with the extent of the distribution because this model predicts that the data of clusters 
should be located only around the point of z co \\ = in Figure 1. 

Finally, we comment on the case of n = —2, which is suggested by the analysis based on the assumption 
that clusters have just formed when they are observed (e.g. Markevitch 1998). Figure 2 is the same as Figure 
lb, but for n = —2. The theoretical predictions are inconsistent with the observational results, because the 
theory predicts rapid evolution of temperatures; there should be few clusters with high-temperature and 
small core radius. 

The results of this Letter suggest fio < 1, and that the clusters of galaxies existing at z ~ include 
those formed at various redshifts. We also show that the model of n ~ — 1 is favorable. Importantly, the 
location of a cluster in these figures tells us its formation redshift. In order to derive the cosmological 
parameters more quantitatively, we should consider merging history of clusters and predict the mass 
function of clusters for each z. 
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Fig. 1. — Theoretical predictions, (a) Radius-density relation (b) radius-temperature relation (c) density- 
temperature relation. Solid line: CIq = 0.2 and M v i r (to) = 10 16 Mq. Dotted line: Slo = 0.2 and 
M vir (i ) = 5 x 1O 14 M . Dashed line: Q = 1-0 and M vir (t ) = 10 16 M Q . Dash-dotted line: Q = 1.0 
and M vir (i ) = 5 x 10 14 Mq. The open diamonds (A, B, a, and b), circles (C and D), and triangles (c and 
d) correspond to z co \\ = 0, z co \ = 1 / O — 1, and z co \ = 0.5, respectively. The observational data (po> R, & n d 
T) are overlaid being shifted moderately in the directions of p and R (/? v i r = 0.06,oo an d r v i r = 8R). 
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